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2. Modification of the large-scale environment by cumulus clouds 

Consider a horizontal area at some le-vl between cloud base and the 
highest cloud top. This horizontal area, shown schematically in Fig. 1, which 
we designate as our unit horizontal area, must be large enough to contain an 
ensemble of cumulus clouds but small enough to cover only a fraction of a large- 
scjle disturbance. The existence of such an area is one of the basic assumptions 
of this paper. 

Because we are not concerned here with acoustic waves, the mass 
continuity equation can be simplified to its quasi -Boussinesq form, 

v ) + &-(/*■) - o. 


where the density p is c function of height only, v is the horizontal velocity, 
V is the horizontal del operator, Vj is the yeriical velocity and Z is the 
vertical coordinate. 

by 

Let &AZ,t) be the fractional area covered the ith cloud, in a 

% A. 

horizontal cross-sec f ion at level Z and time • The vertical mass flux 
through is 


Hj, = J pvs efor — />C5'-I4r £ , ^ 

v, here I cL$ is the area integral over the area <5. and XiK is the average 
vertical velocity of the ith cloud at this level. 

The inward mass flux per unit height, normal to the boundary of the ith 
cloud, is given by from the mass continuity equation (l)(See Fig. 2.) 
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Here the boundary is not necessarily vertical. Then the mass added to the cloud, 

which may be horizontally expanding or shrinking, is 3M-/?Z -f fd&i/d't 

per unit height and unit time. The entrainment and detrainment of mass are given by 


entrainment: B • — ( ■ when y q> • 

*" v o Z dv • , "dZ. ' dt> s 


(3) 


detrainment: 0 L ~ “ when + < q . 




Z> t 


(4) 


E i can be rewritten as )foz + p(d/dt ■k’UT- d/dZ) • Thus, an 

entrainment of mass, which is originally caused by turbulent mixing at the cloud 
boundary, appears either as a vertical divergence of the mass flux within the cloud, 
or as a horizontal expansion of the cloud as it rises, or as a combination of both, 
depending on the dynamics of the cloud. 

The total vertical mass flux by all of the clouds in the ensemble is 

= (5) 

where denotes the summation over all clouds which are p metrating the level 
being considered. 

Let fur be the net vertical mass flux for the large-scale horizontal 
unit area. It satisfies the continuity equation 

vUz) + jz c ^ ) = °j (6) 

where .he bar denotes the average over the unit horizontal area. In general, the 
total vertical mass flux in the clouds, Mq , is not the same as the large-scale 
net vertical mass flux, pur . The difference between Me. and pur is equal 



to the downward mass flux between the clouds. 


M = M c - />w. (7) 

With sufficiently intense cumulus activity, M c can exceed and sub- 

'V 

sidence, O appears in the environment, as in Fig. 1. 

At a given height, some clouds may be detrc.ining while other clouds are 


entraining (see Fig. 1). V/e define total entrainment E and total detrainment D , 
at each level, by 


E- I Ei 

e.c. 


( 8 ) 


and 

0= I o ■ 

d.c. > 


(?) 


respectively. Here 51 denotes the summation over all clouds \ tich are entraining 

at that level, and 51 denotes the summation over all clouds which are detraining 
d.c. 

at that level. E and D , as well as Me / are functions of % . From (3), 

(4), (5), (8) and (9), we obtain 


E- 


D = 


L 

bZ / , 


where 

I 

is the total fractional area covered by all the clouds of the ensemble. 
We define the static energy by 


(10) 


(ID 


s = <v>T + fZ; 


( 12 ) 



where Cp is the specific heat of air under constant pressure , T is the temper- 
ature and Cj is gravity. Cp T is the specific enthalpy of the air/ and gZ is 
its geopotential per unit mass. The static energy $ is approximately conserved 
by the individual air parcel during dry adiabatic processes. Hydrostatic balance, 
which we assume for the environment, gives 



where jo is the pressure , jo 0 is a standard pressure. Ft is the gas constant , 

JL. _ 

and 0 — T is the potential temperature. 

From the budgets of static energy and water vapor in the environment 
we obtain 

^ t [(l-<r e )f s - Es + £0 ;So . 

- U * 

E? + X D-fo. 

where the bar denotes the horizontal average over the large-scale unit area, the 
tilda denotes the value in the environment (which is assumed to be horizontally 
homogeneous for all of the variables) and the suffix Di, the value in the 



detraining air from the ith cloud. L is the latent heat per unit mass of wafer 

vapor. £ is the evaporation of the liquid water detrained from the clouds 

<\/ 

per unit height/ and is ihe radiationa! heating of the environment per 

unit height. is the mixing ratio of v/ater vapor. The first three terms on 
the right of (14) and (15) come from the horizontal area integrations of 
— V'CpV'S) and — over the environment. 

and — represent the convergence of vertical fluxes in the environment. 

Using (10)/ (9)/ (7) and (6) f Eqs. (14) and (15) may be rewritten as 


*\' N ' 

<>S 


(l-«-c)/‘H-£D i (s w - , ?)-L£- Rlr 


d.C, 


“ [v- (pvs) - V-iC'V) s' ] + 


(i -<s-c)f = £ Di (£ w - %-) +■ £ -Myf 

~ [v-Cpxfi - ^*c rz) ?]. 


32 - 


(16) 


(17) 


To simplify Eqs. (16) and (17)/ we use the empirical fact that the total fractional 
area covered by active clouds is small compared to unity. (See 7 for example/ 
Malkus et al./ 1961.) This is consistent with the theoretical finding first obtained 
by Bjerki. (1938), that conditional instability favors the smallest possible 
horizontal cross section for rising motion and the largest possible horizontal cross 
section for sinking motion/ if there is no friction nor entrainment. Asai and 
Kasahara (19 67) found that cumulus convection most efficiently transports heat 



upwards/ and therefore most efficiently releases kinetic energy, when the fractional 
horizontal area of the rising motion is of the orc’er of several per cent. It seems 

that 


cr c< < 1 (18) 

is an acceptable first approximation when we are concerned with only active 
cumulus clouds. 

By definition/ 5 and averaged horizontally over the large-scale 
unit area are given by 


or using ( 11 )/ 


s = c r - <r t j s' + z <Tj s l/ 

(19) 


(20) 

„ /V C" 

s = S + 1 c Si - S> 

(21) 

?= ? + I . 

(22) 


and are the excess static energy and the excess mixing ratio 

of water vapor in the ith cloud over the environmental values. Using the as. jmption 
(18) and the empirical fact that the tempera’ ure difference between the cloud 
and the environment is srnall^ of UWX(Si~ $)« S/ we obtain/ from (21) 


S 


A/ 

S. 


( 23 ) 


However/ the corresponding inequality, £ • c * oes no * 

hold when the environment is not near saturation. Therefore, we use instead 



/ where £*. is the saturation mixing ratio of water vapor in 
the environment. Using this in (22), we obtain 


S' *? ? + <?*-?> «e 

" ?[ l+ J T t<s ‘«]_, P4) 

where V* is the relative humidity of the environment, %■/ 'fc* * When 

)~ CT ’ ' S/ w ^ en the environment is not very dry, the second 

/- _ 

term in the bracket of (24) can be neglected. Then, 

< 25 > 

We further assume that 

V- (/>£) = V- 02 ), (26) 

v- (p\:s) = V- (f £ s ) / (27) 

~^1F2J) ^ V- f'). ( 28 ) 


Here the bar on the right hand sides denotes a running horizontal space average, 
on the scale of the unit area, and not just the average within the fixed area. The 
approximations (26) through (28) are valid when the net lateral horizontal transports 
across the boundary of the fixed large-scale area by cumulus convection (the 
horizontal cumulus eddy transports) are negligible compared to the horizontal 
transports by the large-scale motion. 



Using the approximations (18), (23), and (25) through (28), Eqs. 
und (17) may be rewritten as 



L 

4.C. 


Di(s. r s)-L£-M|| 



is a measure of the static stability of the environment (see (13) and 
(23)), which is usually positive. — represents adiabatic warming of the 

environment when "H -M c “P l ^0(or cooling when M c -pu/< O ) due to the 
vertical motion. For < O (the normal condition) — Pf ^/d% represents 

drying of the environment when M<,-pw>0 (or moistening when M c -pw^<o)^ 
due to the vertical motion. 

In addition to the detrainment and evaporation terms, cumulus clouds 
modify the environment through the r, 'mulus-induced subsidence in the environ- 
ment. The latent heat released within the clouds does not directly warm the 
environment, but if maintains the buoyancy of the clouds against the adiabatic 
cooling due to the upward motion and the cooling produced by the entrainment 
of drier and colder air from the environment. Thus, the latent heat released 
within the clouds main-cans the veiticat mass flux of the clouds and thereby, the 



cumulus-induced subsidence in the environment. The drying nnd warming of the 
environment/ by the cumulus induced subsidence, ore the indirect effects of 
condensation and relca'3 of latent heat, but their vertical distributions can be 
very different from the vertical distribution of the conden$a..on within the clouds. 
This important role of the cumulus induced subsidence in the environment was not 
explicitly made use of in parameterizing cumulus convection until Arakawa (1969)\ 
Eqs. (29) and (30) were derived from budgets for the environment only. 

But they approximately govern the time changes of s and £ , which are averages 
over the total area. This means that the prediction of the large-scale averaged 
field is practically the same as the prediction of the cloud environment, insofar 
as the thermodynamic variables are concerned. This important simp! if i, 
which was used also in the earlier parameterizations by Arakawa (19t>9) and by 
Ooyama (1971), comes from the neglect of accutaf"' e storage of the static energy 
and water vapor in the ensemble of clouds. In fact, we can rederive the right 
hand sides of (29) and (30), as was done by Yanai (1971b, 1973), from budgets 
for the total area. These budgets give 


+ (IQwt«0, < 3 » 

y°ff = -v-(rzp'Jz { rn ) - (£ c r£), 


A description of this parameterization v/as given by Haitiner (1 970, pp. .) 


( 32 ) 



where C- and Q p . are t lie rates of condensation of water vapor and radiational 
heating in the ith cloud per unit height. Using the approximations (27), (28) and 
(26) v/ith Eq. (6), (31) and (32) can be rewritten as 

+ L( 2 C^ 2 ) + ZQ Ri . 




+ © B 


-PZ-vJ , <34 

The quantities inside the brackets are the eddy vertical transports by cumulus 
convection. The eddy transport of S may be written as 

furs -fUrS = ( 21 MiS i + M's) -furS, 

Using (7), (2), (19) and (i 1), the eddy transport can be rewritten as 

puFb -f>&s =Z[ M itSrsJ-GiftiSi-?)] _ 

/V 

The second term in the bracket is negligible compared to the first, when W* '^> WT, 
'here ur is the vertical velocity in the environment, and S -s’. 


Then we have 


fujS -furs ^ Z Ml ($1~ s ) 0 


Similarly, the eddy transport of 9 may be rewritten as 

pH-pH % 2. Mi(gr F>. 

When there is no accumulative storages of mass, static energy and water vapor in 


the cloud ensemble. 




( 37 ) 



( 38 ) 


ES-Z D. st>i - Z MjSi. 

dS. 

^ L XL C L + X Qr i - ° / 

'E ? ~Z D; toi~ fz L %-i (39) 

-ici. =o. 

Using (35) through (39), we can easily show that (33) and (34) c 3 identical to (29) 
and (30), where the heat of condensation and the radiatlona 1 heating in ihe clouds 
do not explicitly appear. 

If. the evaporation of the detrained liquid water takes place at nearly the 
same level where the water is detrained from the clouds. 


t = J- c D t J> Dl , m 

where is- the mixing ratio of liquid water in the air detrained from the i th 

cloud. The above assumption is probably justifiable for the evaporation of 
detrained cloud droplets, but would not hold for rapidly falling raindrops. With 
the assumption (40) and Eq. (7), (29) and (30) can be rewritten as 

/ft = Z C>i(s D1 -L4>i-?)r Me ff 


d, c. > 


-fv-vs -/»«>■ It + e R; 


(41) 


O 


/, H = Z Di($ n + l 0 rf) + M e -| 

■ -/“£•?£ •'/“'ff- . 

The detrainmeni terms will be further simplified in the next section. 


( 42 ) 



3. Budget equations for an individual cloud and assumptions for the detrainment 
We first consider the budgets of mass, static energy, water vapor and 
liquid water for an individual cloud. In the entrainment layer of the i th cloud, 
the budget equations can be written as 


mass: 

M;, ^ 

(43) 

static energy: 

- M -£(M l s t )+LC l + 

(44) 

water vapor: 


(45) 

liquid water: 

M to + Q - R t , 

(46) 


where is the mixing ratio of liquid water in the form of cloud droplets and Rj. 
the rate of conversion of the liquid v/ater to precipitation per unit height. The 
approximations (23) and (25) have already been used. The budget equations in 
the detrainment layer are 


mass: 


(47) 

static energy: 

jzifxsi Si ) = - D £ s w ■ - ) + L C t ■ + Q Rf , 

(48) 

water vapor: 

dt ~ D$oi~ 

(49) 

liquid water: 

0J Di - £(Mi to + Q - R. . 

(50) 


Equations (37), (38) and (39) are obtained by summation of (43), (44) and (45) O'/er 


t 



all entraining clouds and of (47), (48) and (49) over all detraining clouds and 
subsequently dropping the time derivative terms. 

We define the moist static energy by 

5 * = c t>T + + l fr> 

where L ^ is the latent heat per unit mass of the air. /t is approximately 
conserved by the individual air parcel during moist adiabatic processes. We also 
define the saturation moist static energy by 

+ CpT-tJZ-fif* 

where the * denotes the saturation value of the variable. ^0 defines 

the moist adiabatically stable, neutral and unstable lapse rates, respectively. 

We assume that air is saturated in the clouds. Then 


. *5 rcn, p) 1 

- f + df)pCT C -T; 

= r + 

v/here ^ (jTj p ) • Since and k* ~ «S ■+ L g* , 

^ Here, we neglect the effect on of a pressure difference between the cloud 

and the environment. Although the scale analysis . y Ogura and Phillips (1962) did 
not justify *' : s approximation, the recent numerical integration by Wilhelmson and 
Ogura (19/2; indicates that the pressure difference was overestimated by the scale 



where 



For simplicity/ we nea.oct in this paper the difference between temperature and 
virtual temperature/ except for turbulent thermal convection in the subcloud mixed 
layer. Then, the sign of the buoyancy is given by 


buoyancy ^ Q 

according to 


(57) 

Eliminating Cj, from (44) and (45) and from (45) and (46), we obtain 



-Cl _ 

K h l hZ 

+ Qr L , 

(58) 


= E-? - — 
*- or bz 


(59) 


Eqs. (58) and (59) describe the budgets of the moist static e gy and water substance 
for the entrainment layer of the ith cloud. Combining (58) and (59) with (431, we 
mry write 

(&+ w iii)k=~ (60) 

h-D ~ Ri/Mi, ( 6 ,) 

where is the fractional rate of entrainment, defined by 

Mi s ' («) 

The ith cloud may be in its growing stage, with a rising cloud top. We 
then assume that there is entrainment into the cloud at all levels including the 
cloud top. Only after the cloud top has lost positive buoyancy and has stopped 



rising does detrainment take place in a thin layer at the cloud 'op. From (57), 
the level/ 2* at which the cloud top loses positive buoyancy/ is given by 

if 


a t -n 2=li = o. 

Also, from (54) and (55), we have 

( S; - 5 ) 

c 


Eqs. (63), (64) and (65) show that all clouds which lose buoyancy at the same level 
5 share common values, A*(2)j S ( 2) and , for S^and <?. at that 

level . We therefore assume that these clouds are of the same type and have a 

A 

common value £ for Q: at that level. That is. 


A 

is the liquid water mixing ratio at the level of vanishing buoyancy and is not 

necessarily equal to the liquid water mixing ratio of the air which spreads into the 

environment, because an additional condensation (or evaporation) may take place 

a 

nea. .he cloud top due to the concentrated Q R ^there. We let be a function 

A 

of % , but for different levels refer to different types of clouds. 

For the detrainment layer, it is convenient to eliminate Cl f rorn (48) 
and (50). Then we have 

^[/’C l (s l -l&)]= -Di(Soi -Uoi)- 

* L Ri+ Qri . 


Simpson et al. (1965) pointed out that their model supports the rule of thumb 
that the clouds cease vertical grov/th when their temperature soundings recross 
the environment curve. 



We also have the budget equation for water substance/ 

A -£z[v-tt+h)] 

~ Hi . 

Because the thickness of the detrainment layer, / ‘s assumed to be small, 

% 

the mass budget equation for the detrainment layer, (47), may be approximated by 

Dj^Zot = . 

A t- 

Mj at !?r is the mass flux of the. cloud entering the detrainment layer from 
below. Similar simplifications of (65) and (66), and use of (67), (64), (65) and 
(66) give 

C 

— (s - L4 ) + Qri 

and ; t>i.(?t>L + -^0 = D;( •?;)?=£. 

= D c (r+ 2 ) ?= |. . 

The radiaticn term in (68) is retained because the radiation flux at the cloud iop 
* is almost discontinuous, whereas ®g£,^£p£, is finite even when & 2Tp^ 
is infinitesimally small. 

Summations of (70) and (71) over all clouds which are detraining in the 


same thin layer give 


21 D t (s t>l -L^ot)= D(i-l£)+ z Qr, 

^.C, c(.C, 


( 73 ) 



where D is the total detrainment, / _ D . . Substituting (72) and (73) into 

: d.C, 

( 41 ) and (42), we obtain 

= -LDe+M c || +Iq rl 

/'ll = Dtf+£-fr>+ Mef| 

-n-n -r*% %. 

These are the basic equations we use to describe the time changes of the large- 
scale temperature and moisture fields. However, these equations ere valid only 
above the cloud base. The subcloud layer and its interaction with the cloud 
layer will be discussed in a later section. 

Equations similar to (74) and (75), for a three-level model of the large- 
scale temperature and moisture fields, were derived and used by Arakawa (1969). 
Equations which are almost identical to (74) and (75), except for the radiation 
terms/ were derived by Ooyama (1971), Arakawa (1971) and Yanai (1971b). 
Yanai et al. (1973) used these equations to determine the bulk properties of 
tropical cumulus cloud clusters from the observed large-scale budgets of heat and 
moisture. In that study, however, the large-scale tendency terms were obtained 
from observations, while our problem is prognostic and the large-scale tendency 
terms ore exactly what we want to find throjgh the parameterization of cumulus 


convection. 



4. Spectral representation of cumulus ensemble 

Eqs. (74) and (75) clearly show which properties of the cumulus ensemble 
must be parameterized to predict the large-scale temperature and moisture fields. 
The modification of the large-scale fields by cumulus convection depends on: 

(i) the total mass flux in the clouds, Mc(Z) ) 

(ii) the total detrainment from the clouds into the environment, D C30 y 

A 

(iii) the mixing ratio of liquid water at the vanishing buoyancy level, £ (ST,), 

(In addition, cumulus clouds modify the large-scale temperature through their effect 
on radiation.) The problem of parameterization of cumulus convection is now 
reduced to relating these three properties of the cumulus ensemble to the large- 
scale temperature, moisture and velocity fields. 

The total detrainment, D(£^, at different levels refer to different types 
of clouds. When the thickness of the detrainment layer, A'Zpir is infinitesimally 
small, the total detrainment in the layer between 5? and £-4 -Jj? , pCZ)dZ , 

is equal to the total mass flux, at level ^ , of the clouds which lose buoyancy 
within that layer. It is now clear that finding the total detrainment, D(2), as a 
function of height, is equivalent to finding the distribution of the mass flux in 
the different types of clouds which lose buoyancy at the different levels. This 
suggests that we represent the cloud ensemble in spectral form, by dividing the 
ensemble into sub-ensembles, each of which has a characteristic cloud type. 



For simplicity/ we assume that a single positive parameter A can fully 
characterize a cloud type. Then the detrainment level, &p , which is the 
maximum height of the cloud top, becomes a function of A * By choosing X 
properly, we let 5Tp£A) decrease as \ increases, a$ shown schematically 
in Fig. 3. A more specific choice of A will be made later. Let ApC5} 
be X of the clouds which are detraining at level S . is the inverse 

function of ^(A) / which satisfies 

(76) 


identically. Ac>(Z) is also the maximum value of X at level 2 , because 
the clouds which have larger X than A[> have detrainment levels lower than Z . 
The total mass flux in the clouds, M-, can be expressed as 


where 


Me ft)— f JhCZM d\ 

J n ^ 


hvz^)</\= L M c (z) 


(77) 


(78) 


is the sub-ensemble mass flux due to the clouds which have the parameter Aj, in 
the interval (A, 

The total detrainment, OU)d% / ' n fhe layer between £ and £-tc/S$ is 
equal to the sub-ensemble mass flux, at level % , due to the clouds which have 
parameter A^ in the interval AplfcJ-^-^ApfiD/c'Sr) &Z to Apt?) / 
as is shown in Fig. 3. Then we have 




dZ 


( 79 ) 




It is convenient to normalize by 


• h fe,A) & M e i-v ? (za). 

(80) 


(81) 


where is a properly chosen base of the updrafts associated with the clouds. 
Obviously we have 

(82) 

We shall find it convenient to choose the top of the subcloud mixed layer as t ie 
base, , as shown in Fig. 4. Observations/ which we will refer to in the next 
section/ show that the top of the mixed layer is located somewhat below the cloud 
base, z which is approximately the lifting condensation level. The vertical 
mass flux below the cloud base/ then, should be interpreted as the mass flux of 
the updrafts associated with the clouds but not un the clouds. 

Next we consider the budgets of mass, moist static energy and wafer 
substance for the sub-ensemble. The summation of (43) over all members of 


the sub-ensemble and subsequent dropping of the time derivative term give 


ZZ 


2. W) d?. 


Here (78) and (80) have been used. Similarly, (58) and (59) give 


[ •'?(*» Ac 


- ^p- lm. 


5 P 

«[ ^*)(? c «ua>+^< 3!A))]= - s -f rp- _ - r/ (zA} (M) 

Here (78), (80) and (83) have been used. < fr c C£ J ty and are A. , 



sub-ensemble 


*o 

r* 






^ and ■£ , respectively, in the clouds which are members of ihe sub-ensemble. 

YCz *,>) is defined by 

?h(z**)rcz,x )d\~ X Hi , ( 86 ) 

A* e- (AA+tf 

The radiafional heating in the entrainment layer is dropped r «»r simplicity. 

Eqs. (84) and (85) can be rewritten as 

[ hflZM - AtSlJ J < 87 > 

&[it + •&*,*>] = I? c (S?,*; -YlsMSb) 

where ^i|(i? / A)' ls the fractional rate of entrainment for the sub-ensemble, given by 

vkz m 

Except for the cloud microphysical and dynamical processes which determine 
Y'('ZjX) and the subcloud layer processes which determine and .A.) / 

the problem of parameterizing cumulus convec-'on has now been reduced .o ti.iuu.g 
7j , the normalized /ertical profile of the sub-ensemble mass flux, unJ 

, the mass flux distribution frunction at the top of the mixed layer. 

When *f(Z, A) (and, therefore, ) is known, can 

be readily obtained by integrating (84) or (87) with respect to height under given 
!%(%) and ACZqjZS) , and the detrainment level, , and its inverse function, 

Aj>(H,)/ can be found from the condhion of vanishing buoyancy 

h c ^-) ” A. 


(90) 



Cl net 


( s[plZ)\ 

Also, C£,\) . £,.\) can then ^ 3 obta*: d oy infegruiing (85) or (88), wilh 

a parameterized Y'CZ^X) . under given ■> £.(%bA) • j2-C^ls,\) 

is zero because levei is below the cloud bese. From the known 
^CZ + Z\J) anc ^ *' ie vanishing buoyancy condition in the form 


%(Z,AvCi r>)sr 

we can find 

'4(2)= Jtiz,X t> CZ)') y 

which is ihc liquid water mixing ratio at the vanishing buoyancy level. Then, 
only the mass flux distribution function, hzPQ , which is needed fo» computing 
M(2) and 0(2) from (77] and (79), remains unknown* 

Although our knowledge of the dynamics of clouds is far from sufficient, 
the determination of the normalized vertical profile of the sub-ensemble mass flux, 
^ t is logically more straightforward than the determination of the mass 

flux distribution function, . We may assume that the member of a sub- 

ensemble are at random phases in their life cycle and, therefore, the summation 
of the mass flux over all members of the sub-onscmble, as in (78 ) r is proportional 
to the mass flux of a single cloud averaged over its entire lifetime. The constant 
of proportionality i the number of clouds. But the constant of proportionality 
does not matter for ^ (£,A) / since is normalized. A dynamical model 

which governs the life cycle of a single cloud should be able to determine the 



vertical profile of the time-averaged mass flux of that cloud, and, therefore, 

/ for each cloud type characterized by parameter A • However, 
we are assuming that a single scalar parameter A is sufficient to characterize 
the cloud type. For this assumption to be at least approximately valid, we 
must 'hoose A properly. 

The one-dimensional model of the cumulus tower, developed by 
(Simpson et al. (1965) and Simpson and Wiggert (1969, 1971), has been exten- 
sively tested against observations. This model specifies the fractional rate of 
entrainment by 

= A± 

K 

v 

where R is th radius cf the rising cumulus tower and ^ is the entrainment constant 
(see also Simpson, 1971). fR, is either measured or assumed at the cloud base 
and given to the model as an input. The assumption that fc is constant with height, 
in the Lagrangian sense, leads to better agreement with observations than the 
alternative assumptions of horizontally expanding thermals or starting plumes 
(Simpson et al. 1965). 

In our model also, we assume that R* is constant with height in the 
Lagrangian sense. We are not assuming that the cloud is a column-like steady 
jet, and it may consist of a sequence of active elements which have a negligible 
horizontal expansion rate below the level of vanishing buoyancy. But we do 
assume that the fractional raie of enfrainrne rf forth" time averaged mass flux 
of the cloud is approximately •. ans ; lit. We now choose this constant 




fractional rate of entrainment as the parameter X which characterizes the cloud 
type. Alihough the dependence of the entrainment on the radius, as given by 
(94), is not used in this paper explicitly, we may interpret the larger X as 
representing the smaller clouds and the smaller X as representing the larger 
clouds. 


The asumption of constant fractional rate of entrainment greatly simplifies 
the determination of f£,X) . This assumption decouples the determination of 
^2(2 from the so’ution of the entire system cf equations which govern the life 
cycle of a cloud. An assumption about the geometry, like the expanding spherical 
bubbles used by Ooyama (1971), gives a similar simplification. 

Replacing ^UCJz^X) in (89) by X , we obtain 

(95), (82) and the definition of immediately give 


f ~ £p(X) 

0 Zptnik 2 . 

Thus the sub-ensemble vertical mass flux increases exponentially with height due 
to ihe entrainment. Above the detrainment level, CK'iy the mass flux 
becomes zero. Fig. 5 schematically shows for various A * 

To determine ZpC^)^ we must find A C (Z^) • The solution of (84) 
is given by 

2 ? 

A = yj~ [ A A(z> «'] 

Here, (95) has been used. (54) gives 




s c (?,*)-3ar;= v '® 

where is defined by (56). As an examp fe^/icfXA) forgiven h,l'Z) J Jfte, 

and a constani X c CSTg^A) 3 //^ is diown by the broken lines in Fig. 6 and 

<T(*) is shown in Fig. 7. In the figures the pressure coordinate is used. For 

Py-Q (2^0)2 I ar >d ^ ie seconc i term ' n the bracket of (97) vanishes. Then 

^c^O) 5 ® Kjv{ for all % . As A increases/ /t c (£,X) is more reapidly diluted 

by /v(2T) • From (90)/ the detrainment level is given by the intersections of 

the broken lines with the curve /•}*(!?) . The curve /^(A.) , the pressure at 

~th «s obfci Lll-Cd, 

the detrainment level prwliicb. is shown in Fig. 8. The figure shows that smaller 

A 

clouds (larger X ) have lower detrainment levels than larger clouds (smaller A )/ 
because smaller clouds, which have a larger entrainment rate, lose positive 
buoyancy more quickly than larger clouds. 

To find the mixing ratio of liquid wafer at the vanishing buoyancy level, 

A 

which appears in (74), we integrate (85) with respect to height, from a 
given 1 initial ' condition £ c £2? 6/ A) • In addition, a parameterization of the 
rainfall rate, Yt'ZpC) , is necessary. A very crude, but perhaps adequate 
parameterization is being used in the authors' application of this theory (Schubert 
and Arakawa, 1973); Arakawa and Chao, 1973). 

The remaining problems are to find the base level variables !2 e ^ 

C2g,A)f Q nd, most importantly, the mass flux distribution 

function, 7Xo(A). Up to this point our theory is not substantially different from 




h,h*and h c (cal/gm) 






LARGER CLOUDS 


SMALLER CLOUDS 


that of Ooyama (1971) / as far as the basic logic is concerned. He concluded that 
the problem of parameterization of cumulus convection reduces to finding what 
he cal Is / "the dispatcher function" , the rate of generation of buoyant bubbles c;s 
a function of the initial state of the bubbles. However, Ooyama left the deter- 
mination of this dispatcher function an open question and, therefore, his 
paratneterizatic.! was not complete. 


5. Budgets of static energy and moisture for the mixed layer 

* 

In this section we present a model of the subcloud mixed layer, which 
interacts with the cumulus ensemble. Observations over the Carribbean sea by 
Bunker et al. (1949) and by Malkus (1953) show that between the ocean surface 
and the cloud base level there typically exists a mixed layer in which the 
potential temperature, 0 , and the mixing ratio of water vapor, <£. , and 
therefore $ and A. are approximately constant with height. The top of the 
mixed layer is somewhat (about 200 m) lower than the cloud base level. Except 
for the region right below the clouds, there typically exists a very thin transition 
layer, immediately above the mixed layer, in which 0- and, therefore, $ 
rapidly increase and ^ rapidly decreases with height. 

We denote the height of the mixed layer by (see Fig. 9). Zg 
is assumed to be lower than the cloud base, , which is approximately the 
level of lifting condensation. Therefore, we consider only non-saturated mixed 
layers. We model the transirion layer as a discontinuity in S. and Eg • 




In this respect, our approach is similar to those given by Ball (1960), Lilly (1968), 
Deardorff (1972), and Betts (1972). We define 

AS s £ (2e1-) — Sfvjj 

A % - ?(zb + > - ?m, 

-s £ (2 e +) - Am, 

where 5 (Sfe + ), & an< ^ A (J?e + ) are values of S , £ , and A, 

evaluated just above the discontinuity at 2 b and Sh t $ m anc * Am are *' n0 

t 

mixed layer values of S , and A r respectively. aA is given in terms 
o f as and A^_by 

, A A, = A S' A - 

In typical cumulus situations Z.a^_ dominates over as , making aA negative. 
aA* is given by 

. ~ 7Tr7z^) AS ' 

where is defined by (56). Typical vertical profiles of /%, and A* near the 
mixed layer are schematically shown in Fig. (10). Because of the sign of aA , 

' the cumulus cloud updraft prefers to originate from the mixed layer. In this paper, 
we shall only treat that type of situation. 

In the mixed layer 
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where and F^, are the vertical turbulent eddy fluxes of s and ^ . 

Integration of (103) and (104) with respect to Z from zero to Zq gives 

fnYF~ ~ 7 + *e [(h) 0 - [F s ) b ] + (_Qh) M ' C05) 


rM - + t b [( Fj) # -( (>),] 


where 
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(Qr)m = ii I Qz dz, 

J o 

(F s ) 0 and (f^) 0 are the fluxes of s and at the surface, while (F,)^ and 
(Ffi)g are fhe flaxes of 5 and Just beneath Zg (see Fig. 9). The turbulent 
eddy fluxes jump to zero across since fhe turbulence is confined below this 

level. 

To derive equations for the time change of 2g v.v consider the heat 
and moisture budgets in the infinitely thin transition layer shown in Fig. 9. 

Mg is given by 

/ •**>«* 

% e (\) d\ f 


(106) 

(107) 

(103) 

(109) 


(HO) 


which is the total vertical mass flu.: toward the cloud base at level and 
/^/g - (XFg is the subsidence between the clouds at the top of the 
mixed layer. The mass flux into the mixed layer, the thickness of which may 



be 'changing with time, is given by f> B ( D ? * / pt - ) •+ , 

The downward fluxes of S and a* through the top of the transition layer, 
at £=Z e + , are 

(s M + A^) [ •> B (-^ 5 -C?e)+Me] 

L > 


(ill) 


( 112 ) 


The downward fluxes of s and through the bottom of the transition layer, 
at , are 




013) 


£m [^("K- E ~ l,/ e) + M e] ~ ( Vs . 

i 

Here we have defined 

p a . - „ 

Dt ~ XT 

The continuity of heat and moisture fluxes across yields 

' / b -|f= "( M e- / > 8 5 s)- ^?(B) 6 , 

. \ 

This effect is of critical importance when the upper part of the mixed layer is 
saturated and has a layer of stratus cloud (see Lilly (196C)), but not with our 


Old) 


(115) 


(116) 


(117) 



unsaturated mixed layer. Consistency of (116) and (117) requires 

(tye --Z^CMb. 


( 115 ) 


A relation between the fluxes of virtual $ at the bottom and top 
of the mixed layer can be derived from the turbulent energy balance as given 
by Lilly (1968). We lei 

( F sv)e ~ ~ kf FsvJo . 019) 

The flux of virtual 5 is given in terms of the fluxes of S and by 


Fsv = F s + C/>T S Fg. ( 120 ) 

where T is a reference temperature/ wnich is assumed to be constant, and S 
is equal to 0.61. Lilly (1968) discussed the extreme cases k = 0 and 
k “ 1 , which he called the minimum and i. ntrainment cc.es. 

The minimum entrainment case corresponds to total friction dissipation of iu. bulent 
kinetic energy, while the maximum entrainment cere corresponds to zero 
frictional dissipation of turbulent kinetic energy. Dcardorff et al. (1969) 
suggested k=f 0.10, from laboratory experiments on non-steady penetiative 
convection in a water tank. Betts (1972) suggested k £ 0.25. 

(118) and (119) can be combined to give 
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( 121 ) 

( 122 ) 


where 


A5y = AS + Cf>T 


(123) 



Positive ASy is reqjired for dry convective stability. (121) and (122) show 
that with an upward surface flux of virtual £ / there is a downward flux of 5 
and an upward flux of at . This is shown schematically in Fig. 8. 

Eqs. (105), (106), (116), (121), and (122) combine to give 




(124) 

fu !!*=-(« >M • + % [ ( p ll 4 

>0] 

•V 

(125) 

f B ~— = - (Mb' '-e '■'&) + 


(126) 

Eqs. (124) and (125) give 



/„ - (PX )^ <+ f B [ (f a ) 0+ k f£( r s ,)o] 

j> ' (Qr)m^ 

(127) 


where (Fk)o ' s surface flux of K • 

In (124) through (12 7), ASy appears as a denominator. When A5 V is 
small, the direct estimate of ASy may be inaccurate and cause a difficulty. 
Following Deardorff et al. (1969) and Betts (1972), we car. derive the following 
al iernative approximate expression for ASy , which may be used when 
is small (see Appendix I): 

v -k ^ hTi ™ n * s v- < 128 > 

Eq. (126) reflects the fact that what really determines the moss inflow 
into the mixed layer is the entrainment due to turbulent eddies, which depends 
on the turbulent eddy flux of virtual S at the surface. Without the entrainment. 



the top of the mixed layer is simply pushed down by the subsidence Mg'* * 

• . Without cumulus clouds (Mg» O) , the depth of ihe mixed layer 

increases wiiii time, when PgV'/ R + fe (Fsv / )o/‘ 2> S v > 0 • With cumulus cloud.., 
however, the cumulus induced subsidence between the clouds counteracts the 
deepening. When the cumulus ensemble is very active, the subsidence may 
even make the mixed layer shallower. However, the mixed layer cannot be 
made too shallow, beccuse the shallower ihe mixed layer becomes, the less is 
the fraction of air which enters the cloud from the mixed layer and the greater 
is the fraction of air which enters the cloud from the environment above the 
mixed layer; and that environmental air has not been reached by the iurbulent 
upward transport of moisture. A shallow mixed layer is therefore not favorable 
for maintaining an intensely active cumulus ensemble. In this sense, the variable 
%q is a part of the mechanism which controls the total mass flux into the 
- clouds from the mixed layer. Betts (1972) obtained Mg from (126) and (128), 
for a given ur 8 and (F^)© , assuming • Wed o not assume 

^6" Zc 1 however, but we let vary in time, in order to let the 
thickness of the mixed layer be one of the controls on the intensify of the 
cumulus convection. A more quantitative formulation of this mechanism will 
be given later in this paper. 

In this section, we have shown that ^ Sri, %m and, 

therefore, can be determined prognostically. However, whether we can 
represent S ( $c(Zb> A) anc * by the characteristic values 



is a difficult queslion. It depends 


Am 

5 M , and^in the mixed layer, 

on whether the cumulus clouds have their roots in the thermodynamical variables 

within the mixed layer. It has been reported ihat the roots are not observed 

for trade wind cumuli (Bunker et al., 1949; see also Riehl, 1954). The evaporation 

from falling precipitation would give S(Hb,A) lower than S m and 

/ . ^ 

ghcr than , but would leave n^\ unmodified. ■£(%) 

depends on , hut even more on how we parameterize ) r (Z / X~') . 

Apparently, some perturbations of the mixed layer are necessary for triggering 
the onset of a cloud within an otherwise uniform environment. Bui even there, 
the perturbation could be on , rather than on St^ or (see Malkus, 

). We postulate that the primary role of the mixed layer is to supply its 
moisture and static energy to the cumulus clouds, rather than triggering the onset 
of the clouds, and we let 


S c (%8, 

( 129 ) 

%C (%B, ^0 ~ 

( 130 ) 

Ac Am, 

( 131 ) 


However, the existence of more or less organized updrafts at level 
toward the cloud base, is required for the mixed layer to supply its moisture 
to the clouds. 



6. The cloud work function 


Our final problem is to find the mass flux distribution function/ 

we have done up to this point is relatively straightforward. Of 
course, our models of cumulus clouc';., cloud environment, subcloud mixed 
layer, and of their interactions are highly idealized. But oven so, ihe real 
conceptual difficulty in parameterizing cumulus convection starts from this 
point. We must answer the question: how do the large-scale processes 

control the spectral distribution of clouds, in terms of the mass flux distribution 
function, / assuming that they do so at all? This is the essence 

of the parameterization problem. 

In the special case when the mass flux distribution function has a 
sharp maximum oround a certain X , the entrainment relationship (94) will 
give the predominant size of the clouds. In this particular case, therefore, 
finding will also solve the problem of the cumulus size. However, 

as Simpson stated (1971), "Although cumulus size appears to be at least roughly 
proportional to the horizontal convergence in the synoptic regime, what 
really determines the scale of convection remains one of the critical unsolved 
problems in meteorology". 

The solution of }hgO0 may be even more difficult, because we 

must determine the spectral distribution, and not only the predominant size. 

However, in a parameterization theory, it is necessary to find only the most 

probable statistical properties of the cumulus ensemble, under given piacd-' 

Sca/e. coilctCttcsris ^ Q.7\cC 'fip'C tAc £S o-f avi cV? cCiyst <lnaP c/citd 

cut 4. JfLvZ'n fitace. 



and time. Also, what we must obtain is tbe mass flux distribution, and not 
necessarily the population distribution in A"-! 5006 * These two are generally 
not equivalent. 

One might take the view that the mass flux distribution function is 

determined entirely by subcloud layer processes. Such a point of view is 

supported by one-dimensional cumulus cloud models, in which the horizontal 

velocity components and induced pressure gradients are neglected. With the 

exception of precipitation effects, such one-dimensiona! cloud models do 

not allow dynamical interaction betv/een the upper and lower par's of ihe 

cloud. An initial cloud base condition determines the solution only along 

a characteristic line in the £- £ plane. While such a model can predict the 

height and some other properties of the cloud top reasonally well, the prediction 

of the properties of the cloud air which follows the cloud top is doubtful, unless 

proper cloud base conditions are given as a continuous time sequence of 1 initial 1 

1 

conditions. Such ’initial* conditions and, therefore, the time integrated mass 
flux at cloud base can be specified independently from the dynamics of ihe 
cloud, because there is no way for the dynamical processes above cloud base 
to control them, as long as an updraft is produced near the cloud base. Therefore, 
only local subcloud layer processes remain cs a possible mechanism for determining 
the cloud base conditions. 

■The situation is different in models which have more than one dimension. 
An impulse is still needed, if the initial condition is otherwise unform horizontally, 
but only for the first cloud to get started. After the initial time, the cloud base 



conditions for the first cloud and for all subsequent clouds cannot be specified 
but are determined as a part of ihe solution of ihe entire system of equations, 
which includes the dynamics of both the cloud and subcloud layers. Small 
turbulent perturbations in the mix r d layer below the clear environment are not 
likely to trigger new clouds, if ihe top of the mixed layer is sufficiently far 
below the condensation level. Thus, the formation of new clouds between 
existing clouds usually requires stronger impulses, possibly stimulated by ihe 
downdrafts associated with neighboring clouds. Otherwise, a new cloud (or 
a new active part of a cloud) is likely to form in the wake of a preceding 
cloud (or a preceding active part of a cloud), because the solenoida! field 
associated with the preceding cloud (or the preceding active part of the cloud) 
produces a circulation in a vertical plane. 

It was pointed out in section 4 that the sub-ensemble mass flux is the 
population times the mass flux of a single cloud, averaged in time over its 
entire life. At level , we have 

Here )U\)cVl is the population of the sub-ensemble. TdA) is the lifetime 
and Is the vertical mass flux, at level 2 b • of o single cloud 

integrated over its entire life time$. Because all clouds which exist at time t 
must have formed during the time interval ( if "TO') j ~t ) y c/A / T£A) 


(132) 


is the rate of cloud formation. This corresponds to the "dispatcher function" 



of Ooyama (1971)* TUgCA) is ihe foial mass which passes level Zg through 
the entire life of a single cloud. If we represent each cloud by a spherical 
bubble, as Ooyama (1971) did, the toial macs becomes ihc mass of the bubble, 
which depends on its radius only* Because the radius is related to the fractional 
rate of entrainment A through an entrainment relation similar io (94), ihe 
total mass becomes a prescribed function of A , which remains ihe 

same regardless of ihe large-scale conditions. Then ihe mass flux distribution 
function f° r different large-scale conditions is controlled only 

through different "dispatcher functions". This agrees with Ooyama* s conclusion. 

However/ if wc do not assume sphericity, or any other prescribed 
geometry which relates the vertical dimension to the horizontal size, tiro 
functional form for WgfA) is unknown. Wg (?J) is a gross measure 
of the activity of a single cloud of type A , and It is highly probable that 
large-scale conditions control the mass flux distribution function by giving 
different functional forms to r /Y) s (A) • 

The numerical simulation of a cloud by Ogura and Takahaslii (1971), 
with a ‘one and half* dimensional model, clearly shows that the time-integrated 
mass flux near the cloud base is highly sensitive to the rcatc of conversion of 
cloud droplets to rain drops (see Fig. 7 of their paper). When the conversion 
rate is sufficiently small, the cloud attains a steady state ^ (i.c. (X) 

is infinitely large), while the cloud undergoes a life cycle ( is 


1 


The time change of the environment is neglected in this model. 



finite) with larger conversion rates. According to their interpretation, this 
difference is due to the different drag force by rain drops in the middle portion 
of the cloud where the buoyancy force is acting. This does not directly show 
that ‘TVqC/J) is sensitive to the large-scale conditions, but it does indicate 
that £>„) is sensitive to work done by forces in the middle portion of the 
cloud. 


Recently, rapid progress towards the realistic simulation of cumulus 
clouds has been made with two-dimensional models, notably by Arnason et al. 
(1968, 1969), Murray et al. (1970, 1971), Orville et al. (1968, 1970), 

Takeda (1971) and Wilhelmson and Ogura (1972). But so far the studies have 
been limited to the formation and the decay of a single cloud. There has been 
no long-term numerical simulation of a cumulus cloud ensemble analogous to the 
numerical simulation of the general circulation of the atmosphere. 

In this paper, we do not attempt to determine }l( A) , TCX) and (X) 
separately, although that should be an eventual goal of statistical cumulus 
dynamics. But the preceding arguments suggest that v/e look info the generation 
of the kinetic energy of cumulus clouds. 

The time change of the kinetic energy of each sub-ensemble can be 
written as 


dJC(X) 

ett 




( 133 ) 


where and are respectively the kinetic energy and the 

dissipation rate due to all the clouds with fractional entrainment rates between 



A and A+ dX . The first term on the right of (133) is the rate of generation 
of kinetic energy by buoyancy forces. The work done by downward drag due 
to raindrops is included in the dissipation. It is important to note that (133) 
holds for three-dimensional clouds, as far as the clouds are driven mainly by 
buoyancy force. Mechanical interactionswith the verticoi shear of the horizontal 
velocity of the environment and with other types of clouds are neglected. 

A (A) is the kinetic energy generation per unit %(?<■) r/A or the 

efficiency of kinetic energy generation for cloud type A • It is given by 

2p(A) 

AM = J — ^ (s c C2,X)-S(Z>) dz . ( 134 ) 

Be 

.We call A(A) the "cloud work function" for type A • Eqs. (97), (98) and 
(131) give ^ 

(s c a,v-s<z>) = y^- ; [A « + * J ^ztvZewz'- 

for z c CA)< Z £ Zu A), (135) 

cvncL 

,Z 

(5 c (AA}-5^)) «• S M + A | )S(Z') ciz'- 7 (Zj%) SC2) 

J *B 

Joy Z 6 i> Z £ S? t CA).‘ (135) 

The condensation level is given by the consistency of (‘•'35) and (136) at 
% =.2 c CA) » An approximate expression for £ C (A) is given in Appendix II. 

From (134), (135) and (136), we can see that ACA) is a static property of the 
environment, including the subcloud mixed layer. 



In the special case of (A) — £ 3 ^ 


, substitution of (135) 


info (134) gives 



We have the identity 

n M &(*) = a*ov) + / if, //ft'- 1 ) < (i39) 

Using (13?) and (95) in (137,’, we obtain 

- 

-AYVJt / f A (Xfi* 'AV; 

•'Zs* 1 _ .'*e* _ 

f-'Cz) ^(Z) (Am~ h ? *s the buoyancy at cloud base. — bh?'/e'Z 
is positive when the lapse rate is steeper than the moist adiabatic lapse rate. 

When there is no buoyancy at cloud bas° and the environment is saturated 
(A ■= A.*,) f (142) is simply a measure of the conditional instability. 

When the environment is not saturated, the contribution to /\0\) by A — n 
(i.e. L ($-%*")) is always negative and opposes the positive contribution 
by the conditional instability. To have a positive /{(A) , which is necessary 
for kinetic energy generation, the environment must be not only conditionally 
unstable, but also moist enough to give a sufficiently small l A~~ A?' ( 

This effect becomes increasingly important as \ increases. Uncnfraining 


-W'H dz . 


( 140 ) 
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clouds, for which A= O , ore of course not influenced by the environmental 
humidity above the cloud base. 

O can theiefore be considered as a generalized criterion 
for moist convective instability. Because of the entrainment, the criterion 
depends on cloud type. / A) 0 for all A. gives a stable environment. 

7. The quasi-equilibrium assumption 

The purpose of cumulus parameterization is to relate the statistical 
properties of a cumulus cloud ensemble to the large-scale variables, and 
thereby to obtain a closed system of prognostic equations for the large-scale 
variables. However, there is no aprtori reason to believe that this is always 
possible. 

Suppose that initially there are no cumulus clouds but ihe vertical 
distributions of temperature and moisture are such that the cloud work function 
A (A) is positive fora certain range of A • Cumulus clouds will then 
form and develop, and, sooner or later, will enter a non-linear regime. 

In this regime the environment will be modified by the cumulus clouds. 
Typically, the modification will be such that moisture is removed from the 
environment, especially from its lower part, and warming occurs In the 
ens'ironment, especially in its upper part. The cloud work function will then 
decrease in time and eventually approach zero for the entire range of A. 



(i.e. a neutral state is reached). We rail the lime needed for this 
adjur.mcnt to a neutral state "the adjustment time scale", Tadj • After 
this adjusiment is over, there will be no generation of kinetic energy, and 
the clouds will dissipate. We call the time needed for dissipation "the 
dissipative time scales", Tp/s (A) • 

However, if there is c: counteracting destabilization (a generation of 
the cloud work function) by large-scale processes, which we call "ihe large- 
scale forcing", the cumulus activity will be maintained. If the large-scale 
forcing is constant in time, if is probable that the cumulus ensemble will 
approach ari equilibrium. The time needed to reach the equilibrium should 
be about the same as the adjustment time, at least when the large-scale forcing 
is weak. 

When the large-scale forcing is changing in time, the cumulus ensembl 

will not reach an equilibrium. But when the time change of the forcing Is 
s(ourtr~ -than the (icfjust7»et:t time-; 

sufficiently we can expect the cumulus ensemble to approach a sequence 

of quasi-cquilibria. In such a sequence of quasi-equilibria, the large-scale 
forcing and the cumulus ensemble vary in time in a coupled way, and, therefore 
the time scale of the statistical properties of the ensemble is equal to the time 
scale of the large scale fields, • We call this relationship " the quasi- 

equilibrium assumption”. If is also an assumption on paramelerizabilify, 
because, unless a cumulus ensemble is in quasi-equilibrium with the large-scale 
processes, there is no hope that sve can uniquely relate the statistical properties 



of the ensemble 1 o the large-scale variables. 


The quasi-equilibrium assumpiion applied io the specfral disirih.iion 
of the kinetic energy allows us to write 

d KM ^ KC A) 


U t Tls 

From the definition of the dissipative time scale, we have 
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For Tls ?d\$ (*) (tyP' 00 *!)'/ 'O.sMO ^ 5-CC and 

^DtS I/O & I / we see lhat the time change of the kinetic energy 
is negligible in (133), and v/e have a balance between the kinetic energy 
generation and the dissipation. When cumulus clouds of iype A exist 
0h B W > 0) , we have 

cbOO 
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( 141 ) 


(M2) 


(143) 


The right hand side of (143) is the dissipation per unit ))ib 6*0 c ^- 
or the efficiency of dissipation for cloud type . V/e call this the "cloud 
dissipation fund ion". Because the cloud dissipation function must be 
- positive, the cloud work function must also be positive, for any cloud type 
which exists ( 7/7 Q (X) > 0 ). 



8. Generation and destruction of the cloud work function 

In this section, we consider the lime derivative of iho cloud work 
function. Because the cloud work function depends only on the vertical 
distributions of the static energy and the water vapor mixing ratio of the 
environment, including the subcloud mixed la^er, I ho prognostic equations 
*74), (75), (124), (125) and (126), which respectively govern the time 
derivatives of 5 (£) , £ (%) , S M , £ /v} and , are sufficient to jive 
the time derivative of the cloud work function. Those prognostic equations 
involve terms of two types: "cloud terms” which depend on the mass flux 

distribution function either through the total vertical mass flux, Mc(z) (orMp) 
or through lhe total dclrainment, D l£) ) and "large-scale terms" si ch as 
large-scale advection, surface eddy f.uxes and radictional heating leans, 
which do not depend on lhe mass flux distribution !' '.ction. Then the time 
derivative of lhe cloud work function can also be expiessed as a summation 
of cloud terms and large-scale lerms. We may write 


_ fdflW) , B 'X) 

where the suffices c denote the cloud terms. h C\) t which v/e call "la’ge- 
scale forcing", represents large-scale terms. Positive F( A) means genc.ation 
of the cloud v/ork function (destabilization) for type P\, by the large-scale 


(1 


processes. 



The cloud terms (d/\/cit 3c i linearly depend on 
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and D(Z) . Mc(Z> is an integral transform of (A} (see (77) and 

(80). D(£) is also determined by (see (79) end (80)). Thus the 

whole spectrum of cloud types can participate in determining M c (7.) and 
£)(Z) / and/ there r are / in determining ^ c/ A {^O/cLt 3 C . Then (174) 

may be rewritten as 

Aw&X 
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When the kernel $ negative , is the c mount of 

destruction (stabilization) of the cloud work function for typo A through the 
modification of tne environment by typo A * However , the kernel may bocoir 
positive for some combinations of A and X . 

The actual forms of 'P&) and KaX) are given below. A brief 
derivation is given in Appena'ix III. First we define 
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The large-scale forcing is given by 
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where 
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The suffix LS denotes large-scale *cims as previously defined. We cai! 

F C (X) the "cloud layer forcing" and Rv\(X) the "mixed layer forcing". 

The first term in the brocket of (T SO) 7 which is the most dominant term. 


represents generation of the cloud work function due io the cooling cf the environ- 
ment by large-scale processes, typically by adiabatic cooling due to the large- 
scale upword motion. The term which depends on-a/i $~cms to be at least one of 
the dominant terms in (151). When a/l is a large negative, the moist static 
energy in the mixed layer is large compared to that In the environment above the 
mixed layer. The effect of large-scale upward velocity at the top of the mixed 
layer on b'Zftfdt i$ positive, and, therefore, the cloud work function is generated. 
An examination of telafive importance of the terms in the large-scale forcing 
(150) and (151) for actual observed situations will be presented In lh author’s 
subsequent paper (Schubert and Arakav/a, 1973). 

The kernel is given by 


K(A / A.')= K V (AX')+ Ko (A - X )+ KmW), 
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In (154) , the symbol Zf/ has been used for Zp(^) • Ky/ CX ,X) , 
KolW and ktv,(A) originate respectively from M, <^C?) in (74) and 
(75), Df2) in (74) and (75), and K[j in (126). Y/e call Kv(A,A.O the 
"vertical mass flux kernel", k[>(A|X / ) the "detrainment kernel and 
K^CA) the "mixed layer kernel". 

The most dominant term in (153) is the firs! term in the bracket. This 
term represents destruction of the cloud w rk function for type A clouds through 
adiabatic warming of the environment due to the subsidence induced by type 
X clouds. Because f2,A) vanishes above £c>|A) , the upper limit of 
the integral can be replaced by ZplX?) t »f £&(X)> Zpt A7 
Then we can see that Kv (AAO is symmetric with respect to A and 7<f > 



if other terms which follow the first term are neglccied. This symmetry of 
K v i\,x') means that the amount of destruction of A O '* by type Xf 
clouds per unit , (through the mechanism described above) is cqoa 

to the destruction of A^) by type a clouds per unit . This 

situation is illustrated in Fig. 11. 

The detrainment kerne! is always positive for s\ / >s\ • This . ons 
that shallower clouds generate the cloud work function of deeper clouds through 
cooling of the environment due to the evaporation of detrained liquid water and 
moistening of the environment due to ihe detrainment of more moist air from 
clouds* This situation i$ illustrated in Fig. 12. 

The role of varying as a control on ihe intensity of the cumulus 

convection was discussed in section 5. This role is represented by the lest term 
in (1 55 ) f which has a negative contribution to ihe mixed layer kernel because 
is negative. 
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Appendix III 


Derivation of flie Kernels and the Large Scale Forcing 


Using the definition of Aft) given by (134) we can write 
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where ^(7)tx(v) is defined by (146). A term involving ihe lime change 

of 5^(3) does not appear since 2 0 ( A) is a vanishing buoyancy level 

for cloud type X . Using (129), (135), and (136), (lll.l) can be written 
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where ^(T)pix) is also defined by (146). After evaluation of the time 
derivatives in the integrands of (111.2), we can write 
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where f(Z)0- (/!,?-) and ^(4) are given by (147) and (148). The 

f IrsT double integral in (IN. 3) is over the hatched triangular region of 
figure while the second double integral is over the siip'icd area of 'he 


figure. The first double integral can be written as 
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while the second double integral can be written as 
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(III. 5) 


Substitution of (III. 4) and (III. 5) into (III. 3) yields 
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(III. 6) 


Substitution of (74), (75), and (126) into (ill .6) gives 
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Using (77) ano' (79), and changing the detrainmeni term from an integral over Z 
to an integral over X , we obtain 
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(III. 8) 


The vertical mass flux kernel and 1 he mixed layer kernel appear in terms of 
Fredholm type while the detrainment kernel appears in a term of Voltcrra typo* 
- Neglecting the time change of A GO and combining terms we obtain 
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